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ABSTRACT 


This paper considers an alternating renewal process 
with time to failure being distributed exponential (A), 
andawieh asceneral distribution for time to repair. 
Bounds and the precision of the bounds for the avail- 
Bbilvtyerunction were obtained. 

A computer program was written to solve for the 
availability function and some other quantities when the 
down distribution was a gamma (a,8), a integer. The 
availability function with down distribution being a mix- 


ture of gamma-functions is considered. 
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[Peeewecee ETON OF ALTERNATING RENEWAL PROCESS 


A. SYSTEM 

The system considered consists of a single unit or 
Component Which is repaired upon failure and then re- 
turned to operation. 

In the model, the following assumptions are made: 

- Time to failure is a random variable (R.V.) X with 
distribution F (specified below). 

“stime tO periLorm a repair is a R.V. Y with dis- 
tribution G (specified below). 

- Repair commences immediately upon failure. 

- Unit is returned to operating state immediately 
upon compietion of cepair. 

-"Unit wsedsego0d as new after repair. 

- Consecutive operating times between failures (as 
well @s repair times) are independent. 


Tims the system can be modeledeas an alternating re- 


newal process. 


B. DISTRIBUTIONS 
The time to failure - termed uptime - is throughout 
assumed to be exponentially distributed, with probability 


delist tumetion (p.d.f): 


Ge) 


Ae eee) 
Ce). = 0 . tt <0 
) > ee | ee 





Tiree different dastributions for time to repair - 


termed downtime - were treated: 


Gamma distribution. 


Mixture of two exponential distributions. 


- Mixture of gamma and exponential distribution. 


The respective p.d.f.'s are: 


(fe?) 


pol gibt 


g(t) 


(Las) 


g(t) 


@. LAPLACE TRANSFORMS; GENERAL 


Denoting the operating state by 1 and the 


by 0, the objective is to find: 


ct 
N 
° 


failed state 


Sowell: ine “avatlability” of the system, i.e., 


the probability that the system will be "up" or operable 





when called upon at a random time t, given that it starts 
a State L at t = 0. 

Gives toltewing are wadso of interest: 

~aer- the expected number of visits to state 1 
in period (0,t), given that the system is in state 1 at 
c= (). 

- rc probability System iS in state |] at time t, 
Maven 1t Starts in State 0 at t = 0. 

ae seupeceea tine system is "on in’ (0,1). 

it os) shown in [1] that 

(1.3) 
ae) = f Mo tee Xx) 


O 


t 
M, ,(t) il [1 + M,,(t-x)] dG(x). 
O 


Taking Laplace transform on both sides of (1.5) and 
iio, denoting the transform by circumflex ("hat"), i.e., 


M, 5 (s) = f 5 M,; (t)dt, one obtains: 
(1.7) 
ws) = — fees) 
“ s(1-f£(s) -g{s)) 


In similar manner: 
(1.8) . 
P.,(s) Si a 
Sls sires) ) 





(1S) 


Gee ee) eG) 


s(1-£(s)+g(s)) 
ees 0") 


E(s) = ——*(s)___ 
s?(1-#(s)g(s)) 


Only for some of the distributions considered can the 
inversion be accomplished in a closed form. For the gen- 
eral case with gamma down distribution, a computer program 
Pee vileten to pive numerical solutions. 

If a digital computer is not available, good approxi- 
mations are of interest. These - together with bounds - 
have been developed for the availability function. 

Some solutions in closed form in the simpler cases 


are provided. 


—a.© 





II. EXACT SOLUTION - LAPLACE TRANSFORM 


Laplace transforms of the up-and downtime probability 


density functions considered: 




















(202) 

ie Mee. (Gli) 
(2) 

A a 

ets) = aes | MEE (CT 2) 
(5) 

A B, Be 

g(s) =m er + [2-2] B, +s Ref. (1.3) 
(254) 

n ae Bae 2 

g(s) =m Bez + [ten | 2 ss] ; Ref. (1.4) 

1 2 





In the case of (2.2) one obtains upon substituting into 


(i.7) to (1.10) wand rearranging: 





(26S) 
a a 
i, Of 
S*q(s) 
where, 
ea s . 
a(s) = (ste)% +2 Bt (see)? 
i=1 
(7) 
A ae (o} 
BE, i(s) = asus) 
Si gus) 
(2.8) 
a 16 J 
Poi(s) = 
s-q(s) 





(Cure) 


A i: a 
E,,¢s) = 2 2 . 
Ss q(s) 


iieedteto (2.5) and (224), only P,,(s) has been con- 


sidered: 
(2.180) 
m Sa Gas ho» 
Pees) SS ————  ——  ;; Ref. (2.3) 
s[s?+(A+B,+8,)s*B, (A+B) -Am(B,-B,) ] 
C2) 


A (s+B,) (s=B.)- 
i > Ref. (2.4) 
s[s*+As*+Bs+C] 


Piicre 
ee Berge 2B ut A 
Ego, + 22, (A+8,) + ABem 


C = AB*(1-m) + Bae + 248 Bom. 





rt.” INVERSION 


Bic COMMENTS 

Consider the transform of a rational function P(s)/sa(s), 
where the degree of Pie) is less than Salis). This is a suf- 
ficient condition for the use of the residue theorem in 
finding the inverse Laplace transform [2], thus: 


15) 


F(t) = £ (Residues of eSt f(s) at poles of F(s)). 


As is shown in Appendix B only simple poles exist if 
B >. This simplifies Che inversion, 

Equation (2.9) has a double pole at s = Q. The ap- 
propriate adjustment to the procedure used is considered 
leanwe 1 . 

In the sequel, primarily the inversion of P,,(s) will 
be considered. However, since the method used is applicable 
pomcleimcasesewhere Laplace transtorm is of type (3.1), 
the discussion in some detail of inversion of P,i(s) will 


show how the other results were obtained. 


Be INVERSION FORMULA 
Reterrane back to (2.7) 
(forz ) 


Pye) = Lt ese 
S q(s) 


q(s) is a polynomial of degree a. It is shown in Appendix 
A that all roots are complex, except when a is odd, in 


which case one of the roots is rea and negative. 


10 





Paetorine sq(s) into linear terms: 


(5.3). 
Bats) = SCS et) same <3 tee” .) 
where 
we = Be it ee (one= 0 if %o 1s “ould) 
al = Be : Dee jew =i * 5G 


sine twe residue theorem: 








Pee = 2 Residues 

aie lt 

Sie eke 
= lim - 

oie k 

(3.4) 
eS Se: 

B (+8) ” 2s 
2) a ee 

Brad =] pry (ty 8g) 

k=1 

Re= 1.2, ,Omet ek 

Deno t img 


(r,+B)° = Vineet 


o 
math, - te) ges ¥,i 2 # k 
k=1 
then 
(55 ) 
= B ; eee a7 ae Nee 
Py, () ~ Brad : aU Soler byt i Ogre ‘ey? 


Ia) 





(#.oyeont 'd. 


a,t 
ak 


2 2 
a 





sin b, t) 
Ke = 1 oy oe ene = Ls 


Incidentially, for a bounded solution to exist the 
teal part of 211 roots must be negative. 

The numerical solution is contingent on finding the 
roots of sq(s)- Thus can be done algebraically for a < 4. 
The method used is outlined in Appendix C. In addition, 
fom o = 2 the inversion of (3.5) in explicit form is shown in 
Appendix D to illustrate the method. 

For a polynomial of unrestricted degree the computer 
program described in Section V was written. 

If a computer is not available, approximating formula 
Ifepee are prven in Section IV. 

For 1a ae equation (2.9), the double root at s = 0 


leads to: 


2 oe os) 
= Pim eee SS 


s*0 ds s?(s-r,)(s-r,)+++(s-r,) 


- Conse) cr 
ee 
k=1 *°Tk +5 
Th S2S-r 4) 


which upon derivation and letting s > 0 lead to 





(S20) 


Aa(B+%sA(Ca-1)) pone 
(B+aA) * Brad 


E,,(T) = - it 


o} il 


(S20) (EAs) e° 
7) ek 
ie ae Se(s-t Ti(s a ease) cs aa 

hMiscwembar procedwre as outlined for P,,(t) is used for 


the last term (modified by multiplication by s). 


rs 





IV. APPROXIMATIONS AND BOUNDS 
A. GENERAL DOWN DISTRIBUTION 
Consider 


ilegees 


P (SS SS 
a s(1-f(s)+8(s)) 


where uptime X ~ exponential, downtime Y ~ G, unspecified. 








Then, 

__A 

ae A+S 
= 
ial s [ass ats)] 

(4.1) 
Pii(s) = s m 

[ 1-g(s) \ 


By division: 


(4.2) 


P(s) = 


nie 


— E-g(s) , 32 [2-202 ee pso)'s-} 
Ss = =) 


Observing that emis ise is the Laplace transform of the equi- 
librium distribution of G (Reference 4) except for a factor 
Pe |Y)], then 
(4.3) 
P,,(s) = 1 - RG,(s) + R2(G,(s))? - R°(G,(s)) *#°°- 
where 


e BOL; Ges) = HBG, 


tr tm 


14 





With R < 1, convergence is ensured. This implies that 


Dave) > @ et > ~ (Steady state result). Inverting (4.3): 


BGG ies > RG, (t) . p26") (¢) - R*66*) (¢) AGO 
where 6°) (¢) denotes the ith fold convolution of G with 
iecelt. 

heetimg t + ©, 


E[X] 


lim P £)°= ] -9Ree R20eeR? 2. = : 
— 1 () E[XJ+E(Y] 


one obtains the steady state result. 


B. BOUNDS AND PRECISION 

Pemisea reasonable proposition to try to invert the 
first term in (4.3). Inversion of the other terms may be 
chitfreult unless G. is a simple distribution. 

Even though (4.3) is not a series in the time domain, 
Pyeprobabilistic considerations the effect of truncation 
may be deduced. 

Noting that RI g(t) > oe ee (tie then the ap- 
proximation of P,,(t) will give a lower or upper bound on 
Pay it) Wepending on whether the last term is negative or 
Pesttive, €.2.: . 

(4.4) 

ects l= RG, (t) 

(4.5) 

Bait) = 1 - R-Gacey +iReel-? ct). 


ES 





The precision with which the bounds are established 
can be obtained. 

Suppose one trunecateseatter one term. Then the 
maximum value of the remainder is less than R*. Thus, 


(4.6) 
TCE Cc) RGN PE (tL) 21 = Rel (e). 


Peavying found the lower bound for P,,(t), then this value 
mome@ccurate to Within R*. In a rather typical case with 
Pee 1/a0 (ive., steady state = 0.909), the lower bound 
ieee>eablished to within one per cent. 

If it is feasible to compute the upper bound by 
taking two terms in the expansion, then - by similar 
argument - the precision of the upper bound will be R?, 
i.e., to within 1/10 per cent in the case considered. With 
Emretatively large, the precision is poor. If feasible, 


more terms must be computed for a desired precision. 


CG. GAMMA DOWN-DISTRIBUTION 
Consider case with Y ~ Gamma (a,8) a = 1,2,3,°*°: 


(4.7) 





| OL 
nm aed 1 | 2-| (B+s)%-g" 
RG, (s) eee an i: ae = Sap a 
= + 
: s*[3] ee 
Neting that s = 0 is root, one may simplify (4.7) as fol- 


lows: 


(B+s)%-p% _ p%+g% 1o4g% 2¢g+s)5+87 3 (B+5)25+---+(8+8)% 's- 


s?(B+s)° s2(p+s)% 


16 


a 











(4.8) 


a 
f 7a ia 
oo) “23 d F | 





SUMGErtititing imto (4.3): 


(4.9) 


a : 
ra) Ea ae Mio 
at 





Rts 


ae ae} . 


Each term within the bracket contains gamma trans- 
forms, making the inversion simple. 

As an example, consider case with a = 2. Then (4.9) 
becomes: 


(4.10) 

















ee lata} +(e)" ar. 


ite resulting inverse for fixyst approximation: 


(4.11) 


pale oe re, “eee 


B 
Using the correspondence between poisson and gamma dis- 
tributions, the value can be obtained from cumulative pois- 


son tables. 


7 





D. MIXTURE OF GAMMA DOWN-DISTRIBUTION 

Consider a mixture of two gamma-distributions with 
parameters (K,8), (2,y) and weighted (l-m) and m respec- 
mrvety. ihe first term in the expansion peeenes 


(4.12) 


Wir 


P,,(s) ~ 1 - (a H, (s) + 2 f(s) 


where 











A K 
: B Bug)? B 
ics) = ahs -[aSs)* ones (abs 


“ 72 
= it, a eee yi 
HACS) = vag +| she —_ cl 


The first term is seen to be the linear combination of the 


WD 
Ww 








two terms in the down-distribution. For case K = 2, 


R = 3, the inversion is: 


(4.13) 
Beet) > 1. = Mes H,(t) + 24,00] 
where 
E(t) Gu ie ceuliewes (etl) 
H(t) =1-e 8 +1 -eY* (yt+1) +1 - e Y- 


2 

[op sya). 
In case more terms are required, the inversion 

will be somewhat complicated in that there will appear 


cross-terms between the two distributions. The precision 


emcees) will be 


18 





E. WEIBULL-DISTRIBUTION; SePECIAL (CASE 
A lower bound can be obtained when downtime is dis- 
werpucea Weibull (6,2). The reliability functiom is: 
(4.14) 


Wee 
G(t) =e ei A ee ORE 


The expected value is: 


E(t] = ¥= 


B 
tae equilibrium distribution: 


4 2 
ii e Bx dx 


oo 7 > Elta 


which becomes upon re-arranging: 





t 

mM 4X Ag ere i ( oon = (Ze) | 
Yay ~ | cxX ¥ 2 a 

Vor Vi J 

Oo 
(4.15) 
£ 2 

G.(t) “VE J Ae 2) aces 

0 


ines , G, (t) is ~ normal (0,1/28). The lower bound becomes: 
(4.16) 


TT c 1 
ne) > ie ge Ae le) a 


where Z(t/V/28) is normal (0,1). Thus. standard normal 
tables can be used for evaluation. The corresponding 


precision is \*°-1/8. 


19 





i wewieubull-Distrapuevem, Mecreasing Failure Rate 
For the special case that the Weibull distribution 


hes shape - parameter K = 1/p, p = 1,2,3,°°*, the lower 
bound with one term in expansion can be found. The equi- 
Mierium distribution: 


(4.17) 
1 


t = 

7 pP -BzP dz 

G,(t) = pl i e€ 
O 


which upon substitution becomes: 


(4.18) 


i 
lr 
WD 
Le 
' 
WD 
va 
42) 
' 
tf 
Qu. 
~ 


G, (t) 


where 


ae 
Z xX l 


Equation (4.18) is the gamma-distribution in x (=z?) with 


parameters (p,8). The lower bound becomes: 
J 
6S ee ee a CS ED 
gP 
Wath precision: 


Po 
3P 


Ae an example, bet k = %, i.e., p = 2. Thus, 


Pees) > 1 - <3 jeer [ave + 1} | 


The value of term in brackets can be obtained from poisson 


tables. 


20 





V. COMPUTER PROGRAM 


The computer program was written in FORTRAN IV and 
run on the IBM system/360 model 67 computer at the U. S. 
Naval Postgraduate School. The program listing is shown 
in Appendix E. 

The main sequence in the program is: 

Gi) Input data: a,8,A, time increments and number 
of time steps (named respectively M, BETA, UPRATE, TSTEP, 
L in program). 

(ii) Generate coefficients of polynomial in denom- 
Mmereor Of (2.7). 

(iii) Call library subroutine DPOLRT to solve for 
Meets, using Newton-Raphson's method. 

fre} §6<Compute amd store coettirerents for cosine, 
mermccamonin P,,(t), P,,(t)oeMi,(t) and E,,(e); e.¢-, 
see (3.5) and (3.6) 

(v) SummMeqsine coer irerents fer P(t) andwadds ce 
steady state value of P,,(t), to give an indication of ac- 
curacy of final result. Theoretically the sum should add 
Home Fr ilapeing value “TR provided=by™ subroutine to in= 
Gieate aceuracy~of root extraction’should read zero. 

(va) Enter DO-loop for inerementing time and tor 
tonal calculations. 

Gam) Prant, results : t5 Poel te. Po Ct). Sih eee 
Eee, Gt), named : T, PROBUP, DPRBUP, RNWAL, EXPUP respee- 
tively in the program. 


Sampl2 output is shown in Appendix F. 


2 








APPENDIX A: NATURE OF ROOTS 


Polynomial considered: 


(A.1) 
(GIN See eG es RCo Oe) a Sh 
Multiply both sides by x-8, to obtain: 
CA 2) 
Q(x) = xatl - (B-2’)x® ~1p° = P(x) (x=6). 


Let a be even. 
Then by Descartes' Rule of Sign [3], with only one 


variation of sign, there exist at most one real, positive 


root (clearly x = 8). 
Consider 
(A.3) 
Gio = (xl eee ee 


Since there are no variations in Sign, no negative 
real roots exist. CONCLUSION: Fora esr. (A.1) can only 
have complex roots. 

Pemc pe odd. 

By (A 2), one variation indicates the positive Toor 
hoe.  CONSIcer 

(A.4) 


nad a 


@fex) =x ~~ = (BE2) eam ee 
One variation imply one real, negative root. CONCLUSION: 
For a odd, (A.1) has one negative, real root, (a-1) complex 


Oot S.. 


ZZ 





mie roots of (A.1) are ot the form 

Te = a bal ea) ee os lich 2s ae 0 if a is odd. 

Since x = s+8, the roots of s will be 
ges Te Sep 2 aeee baat oof 


and 


a 
se = - X - af. 
Ali’ 


That all a <0 can be veritied by Hurwitz Criterron (3). 
Mewever, since E[X+Y] < ~ and FsG(t) is non-lattice, the 
system will be stable, implying that a; <0 for aly) ay 


(* denotes convolution). 
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APPENDIX B: NON-EXISTENCE OF MULTIPLE ROOTS 


Polynomial considered: 
(Be 1) 


a-1 a-2 


P(x) = x7 +) (x top xX 


ete | 


Ne eibece 


Per simplicity, let a = 6, and assume that a double and a 


eimecle complex root exist: 


bia 


i} 
© 
I+ 


Tz 4 =a fbi 


Ts 6 =) (C> eee 


Consider the usual relationship between roots and coeffi- 


Siemes [3]: 


G5. 2) 
eee oc: ao ; sum one root at a 
time 
(B.S) 
(a+b?) 2(c7+d?) = tee ; sum 6 at a time 
(4.6.5 preduce son 
(Be 4) roots) 
2(a7+b2)+(c7+d2)+4a7+8ac = BA 3; sum roots two at a 
eaies 
Weame Cauchy's inequality on (8.3): 
ae 
(B.5) 
2(a2+b2)+(c2+d7) 3 ie 3 
al ne > (a-+b-)*(c*+d7} Nee 


By (8.4), then 


24 





(B.6) 


2(a2+b2)+c2+d2 = BA - (4a2+8ac) 


Bx 4((ateye-c7) < Bx 


Samee (atc)*-c* > 0. 

Equations (B.S) and (B.6) are in contradiction, since 
3YKB5 > BA for all values of B > X. Thus double roots are 
not possible. The proof can be extended to any a and any 
meltaiplicity of complex roots considered. The proof holds 
also for a odd. 

The two inequalities established for a = 6 hold in 
general ite ett terms of form (a*+b*) will appear in 


isolation from all possible cross terms of form (a°c). 


Zs 





APPENDIX C: POLYNOMIAL DEGREE 3 AND 4 


Provided 8 > X, the roots of polynomial of degree 3 
and 4 can easily be found [5]. 


a= 3 





P(y) = y? + Ay? + ABy + AB? 3 y = S+B. 


Denote: 
a= = (3AB-A2) 
b = 5% (24 °-9A28+27A82) 
1 
A = - 7+ Vb? | am 5 
7°07 
\l 
ee i aes e 
ag || 


Gen the roots, solving for s, are: 





r 
os eb - pp = z (real root) 
ag Wes (A+B) + L (A-B) V3 eee (complex) 
(aos Z Z 5 
a = 4 
Bitxy = x" + Ax? + VRS ee x SB? 3c see 
Denote : 

2 2 
a = = (A-6B8)8 
peau (96-4) Ac6* 
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"A 
m= 2 f 3 


38 = arccos || 
aem 
y =m cos @ + 46/3 
2 2 
R = aa” AB SE Y 
tae 2 2243 ‘ 
D = | 5A? - AB - y + (4A7B8-8AB2-A3)/4R| -. 
nF, 2. 2 
Maen, the four complex roots, solving for s, are: 
oe Se ee 
“| 7 i ii [ia+a] sl di 
ean! Bay ie 1 ‘ 
ie 3 * [i+s| ee ie 


The simplicity in form of polynomial 


Qt] e a 
Piee) = x -~ (B-A}x” - Be 


make solving for roots for a > 4 by manual methods easy, 


for example, by Graeffe's Method [3]. 


fag 





PePENDIX D: AVAILABILITY FUNC€TION IN EXPLICIT FORM; ao = 2 


In order to illustrate the inversion method used, con- 


me Pi, (t) for a = 2: 


Ape) 
sf —- : 2 
Be a ee SS ee 
s((s+8)2+A(s+B) +A8) SiS ris tem) 
where 
r, 5 * 5 Qe oe a 
= + (-A-2B4Ai) - A = YEXB-A2 
Then, 


c 2 SHE 
3 +8)2 e pees O57 ¢ Ses ee 

Eee (Ct) = lin = 3G" se iin BI 
. s+0 Sis T, > Y,J Sean, S\s pI LS 7L od 

it 

(s-r,)(s+B)? en 


lo) Ss aiGae 
S>r, S(S-T, S-T, 


2 2 2 
ase + ty eetal Me Be ete 
i ee rane Gorearer eee sae) 
Since the last two terms are complex conjugate numbers, 
they may be collected and expressed in terms of tricone- 


metric functions: 


B? -@+3)t 
uo) a, 
B7+2X8 


A 
at 
(r,+B)2r,(t,-1,)e- + (eB) rer eee 


YyTo(T,-T2ItT, TT, 


ae 
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B -(B+R)t (actdb)cosst + (ad-bc)sindt 
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